Fingering of Electron Droplets in Nonuniform Magnetic Fields 
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A semiclassical analysis of a two-dimensional electron droplet in a high, nonuniform magnetic field 
predicts that the droplet will form "fingered" patterns upon increasing the number of electrons. We 
construct explicit examples of these patterns using methods first developed for the flow of two- 
dimensional viscous fluids. We complement our analytical results with Monte Carlo simulations of 
the droplet wavefunction, and find that at the point where the semiclassical analysis predicts a cusp 
on the interface, the droplet fissions — a type of "quantum breakup" phenomenon. 

PACS numbers: 73.43.-f,05.10.Ln,47.54.+r 
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The growth of many physical systems is limited by 
diffusion — e.g., the diffusion of heat for solids growing 
into supercooled liquids hJ, the diffusion-limited aggre- 
gation (DLA) of colloids |3Lor the diffusion of magnetic 
flux in a superconductor y. The unstable growth of 
all of these systems Q produces beautiful patterns. An 
important paradigm for this class of pattern- forming sys- 
tems is Laplacian growth (LG) — diffusion-limited growth 
in the limit of long diffusion lengths Agam et al. 

01 recently discovered an intriguing connection be- 
tween LG and the growth of a two-dimensional (2D) elec- 
tron droplet in a high, nonuniform, magnetic field. Their 
semiclassical analysis of the droplet wavefunction reveals 
that as the number of electrons in the droplet is increased 
the droplet maintains a uniform density while its bound- 
ary evolves according to the LG equations. The unstable 
growth leads to a "fingered" pattern! 

We expand upon the work of Agam et al. in several 
important ways. First, we show that a simple magnetic 
field inhomogeneity outside the growing droplet yields 
tractable analytical results for the droplet growth which 
allow us to predict the onset and structure of interfa- 
cial singularities. Second, using an analogy between the 
probability density for the electrons and the Boltzmann 
weight for a fictitious 2D plasma in a background po- 
tential provided by the magnetic field inhomogeneity, 
we perform Monte Carlo (MC) simulations iTj] of the 
droplet growth. The droplet breaks apart beyond the 
point where the semiclassical analysis predicts a singular- 
ity, providing a quantum counterpart for droplet breakup 
in classical fluids Finally, using the MC method we 
follow the evolution of the droplet interface in a random 
magnetic field and find an intricately fingered pattern. 

Hele-Shaw flow provides a simple example of LG. Two 
glass plates confine air and water, with the air injected 
into the center at a constant rate and the displaced water 
extracted at the edge, see Fig. The air domain is at 
a constant pressure p = and the velocity of the water is 
determined by Darcy's law, v = — Vp; since V-v = 0, the 
pressure in the water is a harmonic function of z = x + iy 
with a sink at infinity. On the interface w„ = —dnP, 
and in the idealized version of LG the surface tension is 




FIG. 1: Growth of a 2D electron droplet in an nonuniform 
magnetic field; on the left the interface is shown for succes- 
sively larger electron number. In the Hele-Shaw analogy, the 
electron droplet corresponds to air which displaces water. 



neglected and p = on the interface. Conformal mapping 
methods [sl, may be used to map the exterior of the 
unit circle in an auxiliary w-plane onto the water domain 
in the z-plane at each instant of time t using z = z(w, t) 
(see Fig.0, with IleiwdwZ dtz) = 1/2 on the interface. 

Next, consider a flat interface and produce a sharp 
outward bump of the air into the water which locally 
compresses the water isobars and produces a large pres- 
sure gradient near the tip. Then w„ = —dnP implies that 
outward bumps grow more rapidly [l]| than the flat in- 
terface, and absent a stabilizing effect such as surface 
tension this bump develops into a singularity in finite 
time 0,01 for most initial conditions. 

We next turn to the many-body wavefunction for 2D 
electrons in a perpendicular, nonuniform, magnetic field. 
If we ignore Coulomb interactions and scale magnetic 
fields by the average field Bq, lengths by Ih = ^Jhc/eB^, 
and energies by heBQ/2mc, then the Hamiltonian is 



H 



n- 

i=i 



-iVj+A,)' + ^a,B{r,) 



(1) 



If g — 2 then the many-body ground state wavefunction 
for Eq. can be found exactly even for a nonuniform 
magnetic field [l^ ll^[l4 ; for a spin-polarized, filled Lan- 
dau level — 1) |5l.ll3l 1^ we have 
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where W{z) = -|zp/4 + V(z), with V{z) a "potential" 
which solves V^V = -SB = -[B{z) - BoI/Bq, and tn 
is a normalization integral. Following Ref. |5( we assume 
that SB = in the region of the droplet; V{z), however, 
is nonzero in this region and the electrons are influenced 
by distant field inhomogeneities — a manifestation of the 
Aharonov-Bohm effect In the region of the droplet 
V{z) is an analytic function and has the expansion in 
terms of the harmonic moments tk [l^ 



F(z) = iRe^^,z^ 
fc=i 



tk 
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Although the wavefunction (O is exact for the nonin- 
teracting system, extracting physical observables such as 
the density is a formidable task. Agam et al. 5] perform 
a semiclassical analysis of tm for ^ 1 and find that 
the integral is dominated by configurations in which the 
electrons are uniformly distributed in a domain of area 
A — 2ttN with moments tk given by Eq. 0. Further 
analysis establishes the connection with the LG problem 
and provides us with powerful techniques (lOj for deter- 
mining the density distribution. 

The derivation of the wavefunction ISJ requires two im- 
portant assumptions — noninteracting electrons and g ~ 
2 — which deserve a brief comment. First, for a uniform 
magnetic field Laughlin's incompressible liquid state does 
an excellent job of capturing correlation effects due to 
the Coulomb interaction 0. However, this may not be 
correct when the magnetic field is nonuniform as the 
Coulomb interaction smoothes density fluctuations and 
will reduce the fingering of the droplet, tending to make 
it circular. In fact, one can think of the Coulomb in- 
teraction as introducing an effective surface tension into 
the problem 0|. In this case one can speculate that 
the fingering only occurs as a driven, nonequilibrium 
phenomenon — the droplet must be constantly "pumped" 
by adding electrons to cause the droplet to finger and 
evolve [l7|. The dynamical effects of Coulomb interac- 
tions for a classical 2D electron droplet were studied in 
Ref. perhaps this approach could be merged with 
the work of Agam et al. in the semiclassical limit. The 
Coulomb interaction could also cause the droplet edges 
to reconstruct jl9j . which may also alter the fingering. 
Second, g 7^ 2 in semiconductor devices, and one might 
worry about the spin polarization of the droplet. In uni- 
form fields the Coulomb exchange interaction stabilizes 
the polarized 1^ = 1 state even for g = [l^, and it is 
likely that the same stability occurs in nonuniform fields, 
which would allow for a perturbative analysis in g — 2 
[ll^ . It is unlikely that g 2 will produce any qualita- 
tive changes of the wavefunction ^ and its analysis. 

As a concrete example we will consider field inhomo- 
geneities which produce a monomial potential of the form 



V{z)^ -ReitM+iz^'+'). 
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with M a positive integer. This is the leading term of 
an interior multipole expansion of SB. The beauty of 
this choice for V{z) is that the corresponding conformal 
map which describes the evolution of the interface in the 
semiclassical limit has the simple form 
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where r and um are real functions of N . Using either 
the moment method described abov e l21l o r the Schwarz 
fimction method discussed in Refs. [lO|,l23i ^^e parame- 
ters in Eqs. Q and (0) can be related as 
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t = 2N 



MuIj. 



(6) 



(M + l)r^-f 

Maps of the form ISJ have been studied ^3 for the Hele- 
Shaw problem, and many of the results can be directly 
applied to the electron droplet. For M = 1, the interface 
is an ellipse which evolves smoothly, with constant eccen- 
tricity, for increasing N [2^. The behavior for M > 2 is 
more interesting; the harmonic measure w'{z) may have 
poles which coincide with the interface at some critical 
value of the electron number N*, resulting in singular- 
ities. A detailed analysis yields the following: (i) The 
interface develops M + 1 simultaneous cusps at 

N* = i(M-l)M-(^-^+i)/(*^-i)[(M + l)iM+i]"'/(*'"''. 

^ (7) 
(ii) If X* = x{N*) is the position of the cusp and x{N) 
is the position of the finger for N < N* , then 



x{N) = 



y/N* - N + 0{N* - N), (8) 



so that the approach to the critical value is univer- 
sal, with a "velocity" Vn = dx/dN which diverges as 
{N* — A)^^/^. (iii) If we define the scaled variables 
X = {x* — x)/x* and Y = y/x* , then the shap e of the 
cusp is universal and has the "3/2" form 
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These interfacial cusps arise from the neglect of surface 
tension in the idealized LG problem. For the Hele-Shaw 
problem the surface tension smoothes the cusp, while for 
the electron droplet one expects that fluctuations beyond 
the semiclassical approximation will smooth the droplet 
Q. However, the mere existence of the singularity sug- 
gests that something dramatic will happen to the droplet, 
and to investigate this we need to go beyond the semi- 
classical level and obtain an essentially exact evaluation 
of the particle density using the Monte Carlo method. 

The probability density for the electron droplet can be 
written as \^\'^ = exp{-pUci)/NlTN , with 
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which is the Boltzmann weight for a classical 2D plasma 
in a background potential —V{z), with a partition func- 
tion Tpf and P = 2. The dominant contributions to r^r 
can be determined using the Metropolis algorithm, anal- 
ogous to Laughlin's studies 7] of the quantum Hall effect 
in a uniform field [i.e., V{z) — 0]. An initially random 
distribution of N particles is equilibrated by moving each 
particle a distance of 2.0 in a random direction, with 
moves accepted or rejected according to the detailed bal- 
ance condition; the average acceptance rate for the moves 
is close to 0.5. After an equilibration period of about 10'^ 
moves per particle we compile a histogram of particle po- 
sitions in bins of size 0.3 x 0.3 for about 4 x 10'' moves per 
particle, and then calculate the average density in each 
bin. The electrons were confined by the walls of an im- 
penetrable "simulation box" whose dimensions were cho- 
sen to be large compared to the linear dimension of the 
droplet. We accurately reproduce known results for 
the droplet size, density, and energy for a uniform field. 
For a nonuniform field the density forms a complicated 
pattern, but it still has a uniform value pavc well away 
from the interface. We determine the interface by finding 
the locus of points whose density is between, say, 0.3/5avo 
and 0.4pavc (other definitions give similar results). 

We have performed numerous MC simulations for or- 
dered and random magnetic field configurations, with the 
results shown in Figs. Q-Q. Let's start with the poten- 
tial V{z) = (^3/2)2^ for M = 2, which may be thought 
of as the leading term in a multipole expansion for an 
arrangement of six thin solenoids with alternating flux 
±$ placed on the vertices of a hexagon a large distance 
R from the center of the droplet, so that — 4$/i?'^. 
The semiclassical analysis predicts that three cusps will 
appear at the critical electron number N* — l/(144i|), 
and we have chosen N* = 500 for these simulations. The 
excellent agreement between the semiclassical and MC 
results can be seen in Fig. |(5J); for TV = 494 the nascent 
MC cusp is slightly rounded compared to the semiclassi- 
cal prediction due to the discrete nature of the particles. 
We find equally impressive agreement with the position 
of one of the tips, Eq. ((SJ, and the shape of the cusp, 
Eq. jnj. When N > N* the electrons in excess of N* 
move toward the edges of the simulation box as shown 
in Fig. Q, so that the central droplet fissions. In the 
2D plasma language, the incompressible fluid of particles 
saturates the potential minimum and additional parti- 
cles which are "poured" into the potential "spill out" 
and accumulate at the edge of the simulation box. Fi- 
nally, Fig. Q shows a simulation for a random distri- 
bution of thin solenoids, with a highly ramified pattern 
for the averaged electron density. Determining whether 
such a pattern is fractal requires a study of the scaling 
law N ~ R^f , with R the radius of gyration and D f the 
fractal dimension (with Df — 2 for compact and Df < 2 
for fractal patterns). An unambiguous determination of 
Df requires simulations over many decades in N; this 
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FIG. 2: Comparison of semiclassical (red) and Monte Carlo 
results (blue) for the droplet interface with M = 2; starting 
from the innermost curve, iV = 50, 100, 200, 300, and 494. 
Three cusps are predicted to form at A'^* = 500. 




FIG. 3: Snapshot of particle positions in a simulation for 
M = 2, A* = 500, and A = 600. Eighty nine electrons have 
split from the central droplet and are organized into three 
smaller droplets at the boundary of the simulation box. 

is currently beyond our computational resources and we 
leave its resolution as an interesting open problem. 

In summary, we have constructed several explicit ex- 
amples of magnetic field inhomogeneities which cause a 
noninteracting 2D electron droplet to finger. The semi- 
classical results (idealized LG) are in excellent agreement 
with our Monte Carlo simulations of the droplet wave- 
function. An important implication of this result is that 
if there is an effective "surface tension" in the nonin- 
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FIG. 4: Average electron density for 280 electrons in a ran- 
dom array of 400 solenoids; blue is low density and red is 
high density. In this simulation 50 "attractive" solenoids were 
placed far from the droplet, with 350 "repulsive" solenoids 
placed closer to the droplet. The total solenoid flux is zero. 

teracting electron droplet problem it must be extremely 
small, and the regularization of the cusps appears to be 
controlled by the discreteness of the particles. Beyond 
the semiclassical singularity we find that the droplet fis- 
sions. Recent advances in imaging techniques for the 2D 
electron gas Q may make it possible to observe these 
fingered structures; an easily tuned field inhomogeneity 
could be provided by placing the electron gas in close 
proximity to a type-II superconductor in the vortex state 
[l3j. Our results also raise the interesting question of 
whether equilibrium Monte Carlo methods using a ran- 
dom V{z) can be used to efficiently simulate nonequilib- 
rium growth processes such as DLA 

3; 

in other words, 
what can the quantum Hall effect tell us about DLA? 
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